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ABSTRACT 


This  report  presents  a  preliminary  study  of  the  sensitivity 
analysis  for  dynamic  systems  with  emphasis  on  its  applications  to 
structural  control.  Definitions  are  first  given  for  different 
sensitivity  functions  in  the  time  and  the  frequency  domains.  Since 
most  physical  quantities  of  dynamic  systems  cannot  be  expressed  in 
analytical  forms,  we  introduce  an  indirect  approach  to  determine 
their  sensitivity  derivatives  from  the  sensitivity  equations 
derived  from  governing  equations.  A  direct  application  of  the 
sensitivity  analysis  can  be  found  in  the  integrated  control  and 
optimization  in  which  design  variables  and  control  variables  are 
treated  equally  as  the  system  parameters  active  in  optimization. 
An  extensive  review  and  evaluation  of  the  existing  techniques  in 
this  area  are  given  to  identify  a  feasible  algorithm  for  future 
improvements.  Finally,  a  new  control  algorithm,  called 

optimization  based  instant  control,  is  proposed  for  those  systems 
subjected  to  general  deterministic  or  random  excitations.  Unlike 
the  conventional  algorithm,  the  optimal  control  is  designed  and 
implemented  according  to  instant  information  of  the  excitations. 
The  important  feature  of  this  approach  is  that  the  original  optimal 
control  becomes  a  problem  of  static  parameter  optimization.  The 
formulation  layout  makes  it  possible  to  apply  the  newly  developed 
compound  scaling  algorithm  [20]  in  optimal  structural  control. 
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1.  INTRODUCTION 


Sensitivity,  as  the  term  suggested,  is  of  concern  with  the 
effects  of  the  variation  of  the  parameters  to  the  system  response 
and  characteristics.  This  concept  was  first  introduced  by 
mathematicians  to  evaluate  the  behavior  of  the  solutions  for 
differential  equations  affected  by  their  coefficients  [17].  Early 
application  of  this  technique  was  mainly  in  the  assessment  of  the 
discrepancy  between  the  actual  system  and  its  theoretical  model. 
The  situation  has  drastically  changed  during  last  three  decades. 
Growing  interest  of  the  sensitivity  analysis  has  widely  expanded  to 
many  research  areas.  In  particular,  many  sophisticated  and 
optimization  theories  have  emerged  recently  from  this  important 
technique. 

The  optimal  design  is  the  approach  of  searching  a  set  of 
system  parameters  that  correspond  to  minimum  values  of  certain 
objective  functions.  The  parameter  search  is  guided  by  the 
sensitivity  derivatives  of  those  functions,  which  determine  the 
steepest  gradient  directions.  How  to  determine  the  sensitivity 
derivatives,  hence,  is  essential  for  the  optimization.  For  the 
case  of  static  loading,  theory  has  been  well  developed  and  research 
efforts  have  been  concentrated  on  the  development  of  more  efficient 
numerical  techniques  [18,19,1],  For  the  case  of  dynamic  loading, 
complication  arises  in  the  parameter  optimization.  This  is  due  to 
the  fact  that,  in  general,  a  structural  response  cannot  be 
expressed  in  a  closed  form  and  many  physical  variables  are  time 
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dependent.  It  is  conceivable  that  the  determination  of  the 
sensitivity  derivatives  for  such  a  case,  if  possible,  may  lead  to 
a  very  cumbersome  process.  Thus,  the  possibility  of  calculating 
the  sensitivity  derivatives  and  the  efficiency  of  numerical 
computations  become  two  important  issues  for  the  researchers  in  the 
area. 

Application  of  the  sensitivity  analysis  in  the  control  was 
first  performed  by  Bode  [2] ,  who  established  different  sensitivity 
measures  in  the  frequency  domain.  These  sensitivity  measures 
became  very  important  in  the  classical  control  theory.  However, 
the  classical  control  theory  has  limited  applications  for 
structural  system  since  most  of  the  structures  have  multiple  inputs 
and  multiple  outputs.  Structural  control,  thus,  must  be  pursued 
using  modern  control  theory  or  optimal  control  theory.  The 
sensitivity  analysis  again  is  playing  an  increasingly  important 
role  in  the  modern  control  theory  [5]. 

Recently,  Hale,  et  al.  [7]  and  Haftka,  et  al.  [6]  discussed 
the  possibility  of  improving  control  performance  by  choosing  proper 
system  parameters.  They  suggested  that  the  control  and  parameter 
optimization  be  performed  simultaneously  in  order  to  achieve  the 
best  control  design.  Many  research  works  along  this  line  have 
shown  promise  results  for  the  case  where  the  structure  is 
experience  zero  or  white  noise  excitations  [8,9,13,14].  Very  few 
have  shown  their  successes  for  cases  otherwise.  One  of  the 
objectives  of  this  preliminary  research  is  to  evaluate  existing 
techniques  in  the  integrated  control  and  optimization,  and  identify 
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the  most  feasible  approach  to  be  adopted  in  the  future. 

It  is  known  that  the  optimal  control  theory  is  established  on 
the  basis  of  optimization,  in  which  the  sensitivity  derivatives 
play  very  important  roles.  However,  the  conventional  optimization 
scheme  is  not  necessary  if  the  external  excitations  are  zeros  or 
white  noises.  It  can  be  proved  that  the  optimal  design  for  such 
cases  are  corresponding  to  the  solution  of  the  Riccati  equation,  a 
nonlinear  ordinary  differential  equation  with  closed-form 
coefficients  [11].  In  reality,  however,  most  external  excitations 
applied  to  structures  are  not  white-noise;  thus,  the  conventional 
optimization  scheme  must  be  adopted.  In  this  report,  we  propose  a 
new  approach,  called  optimization  based  instant  control,  which  is 
valid  for  structures  subjected  to  general  excitations.  The  main 
feature  of  this  approach  is  that  one  only  needs  to  solve  a  static 
parameter  optimization  problem  since  the  optimal  control  problem 
can  be  converted  into  this  form. 
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2.  SENSITIVITY  ANALYSIS 


2.1  Definitions 

Consider  a  real  or  complex  function  F(a)  with  parameter  a. 
The  sensitivity  function  at  the  nominal  value  a  =  a0  is  defined  as 


5_  ma) 

da 


(2-1) 


which  indicates  the  rate  of  change  of  the  function  F  at  the  nominal 
value  a0.  The  reason  of  using  partial  derivative  is  that,  in 
general,  F  also  may  be  function  of  other  parameters,  or  variables. 
This  definition  only  refers  to  the  sensitivity  with  respect  to  a. 
Equation  (2-1)  is  also  called  an  absolute  sensitivity  function. 

Now  we  extend  the  definition  for  the  sensitivity  to  the  case 
of  r  parameters,  indicated  by  the  vector  a  =  [at,  ...  ar]T.  The 
sensitivity  of  F(a)  is  characterized' by  a  sensitivity  vector  S  with 
the  j-th  component 


C  _  dF(  a) 

r^7 


j  =  1,  ....  r 


(2-2) 


Apparently,  equation  (2-2)  is  identical  to  equation (2-1)  except  for 
a  being  a  vector. 

In  a  general  case  where  the  system  response  or  other 
characteristics  are  described  by  a  vector,  say  F(o)  =  {Fj(a)  ,  ~, 
Fn(a)}T,  we  define  an  n  X  r  sensitivity  matrix  S.  of  which 
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<H>) 


3a,  l*"‘ 


j  =  1,  •••,  r;  k  =  1,  •••,  » 


(2-3) 


Equation  (2-3)  represents  the  sensitivity  of  the  k-th  components  of 
P  with  respect  to  the  j-th  component  of  parameter  vector  a. 

In  the  classic  control  theory,  it  is  useful  to  employ  the 
relative  sensitivity  function  in  the  frequency  domain  [2].  The 
relative  sensitivity  function  is  defined  as 


Q  _  d  InF(a)  ■ 

d  lna  lo=# 


(2-4) 


The  relativity  of  this  sensitivity  function  may  not  be  obvious, 
unless  we  rewrite  equation  (2-4)  as 


_  3F(a)/F( a) , 

0  ‘  3«/a  I™* 


(2-5) 


which  represents  the  ratio  of  the  percentage  change  of  function  F 
to  the  percentage  change  of  parameter  a.  The  relative  sensitivity 
function  is  related  to  the  absolute  sensitivity  function  by 


Q  =  S 


^(«o> 


(2-6) 


If  P  and  a  are  both  vectors,  the  relative  sensitivity  is  defined  by 
a  rectangular  matrix  Q  with  elements 


o  -  c  F*(a)i 


j=l,-,r;k=l,-,n 


(2-7) 
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2.2  General  Rules 


It  has  been  seen  that  a  sensitivity  functions  are  the  partial 
derivatives  with  respect  to  a  parameter  vector.  The  general  rules 
thus  are  basically  the  rules  for  the  partial  derivatives.  For  the 
sake  of  the  completeness,  we  list  these  rules  in  the  following: 
a .  summation  rule 


Sa  =  sal  +  sa2 


b .  product  rule 


F  F  F  F 

s*  2  -  VW  + 


c.  chain  rule 


“a  - 


d.  quotient  rule 


SFi  5«‘  Fl(a0>5« 


F2(ao)  F22(ao) 


(2-8) 


(2-9) 


(2-10) 


(2-11) 


where  a0  is  the  nominal  value.  In  the  above  equations,  the  super- 
index  indicates  the  function  for  which  the  sensitivity  is 
calculated,  and  the  sub-index  indicates  the  parameter (s)  with 
respect  to  which  the  sensitivity  function  is  measured.  For  the 
relative  sensitivity  function,  we  also  have 
e .  product  rule 


f.  quotient  rule 


g.  chain  rule 


(2-13) 


(2-14) 


2 . 3  Output  Sensitivity 

It  is  important  to  note  that  the  sensitivity  functions  can  be 
determined  directly  from  their  definitions  provided  in  the  last 
section.  This  is  not  the  case  for  dynamic  systems  in  which  the 
objective  functions,  such  as  state  variables,  are  involved  in 
differential  equations.  The  objective  function  is  implicitly 
given.  Therefore,  we  must  look  for  other  alternatives  and 
determine  the  sensitivity  functions • indirectly. 

Consider  a  system  governed  by  the  differential  equation 

F(y(H\-,y(1)y,u,aj)  =  0;  £_15) 

y°\0)=yQ\  j=l,-,n 

in  which  F  is  a  nonlinear  function,  y(t)  an  output  function  and 
u (t)  an  input  function.  The  numbers  within  the  parentheses  indicate 
the  order  of  the  derivatives  respect  to  time  t.  Our  goal  is  to 
calculate  the  sensitivity  functions  S£,  j=0,l,...,n,  at  the  nominal 
values  a?,  j=0,l,...n. 

Taking  derivatives  with  respect  to  otj  on  both  sides  of 
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differential  equation  (2-15)  results  in 


BF  By  (B>  ( 

-  dF  ,SF 

Byin)  d<ij 

QyW  da.  daj  3a j 

**|  =o- 

Ba.  I,=0  U’ 

** 

II 

H 

\s 

* 

<*  • 

II 

(2-16) 


Use  is  made  of  the  fact  that  input  u  and  initial  condition  are 
independent  of  parameter  vector  o.  Recognizing  that  derivatives 
with  respect  to  t  and  atj  are  exchangeable,  we  evaluate  equation  (2- 
16)  at  the  nominal  value  a0,  and  obtain 


in  which 


(2-17) 


dyfog) 

da . 


(2-18) 


It  is  interesting  to  see  that  the  sensitivity  equations  represent 
a  time  variant  linear  system  in  spite  of  the  original  governing 
equations  being  nonlinear. 

For  a  special  case  where  the  system  is  governed  by  a  linear 
time-invariant  differential  equation 
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bH(a)y(tt)  +  -  +  bx(  a)ym  +  b0(a)y  =  u(t); 

y(k)(0)=y£);  k=l 

The  corresponding  sensitivity  equation  is 


(2-19) 


KS?  ♦  -  ♦  \sf>  |.o'<1V>) 


sf\ 0)  =  0 


k=l ,-,n;  j=l,~,r 


where  y(a°,t)  is  the  nominal  solution  of  equation  (2-19), 


y(a  V)  =  f  f/i(aV-t)w(T>iT  (2-21) 

JO 

in  which  h(a°,t)  is  an  impulse  response  function.  Equation  (2-20) 
is  identical  to  the  governing  equation  (2-19)  except  for  the  input 
term  on  the  right  hand  side.  Solution  for  equation  (2-20)  can  be 
written  as  a  convolution  integral;  namely, 


sr  /;*(«“, T-t)[^  |„„yW(T,«») . ...  ♦  |l  |„.y«>(a0,T) 


(2-22) 


dbn 

+  -r-2-L»y(a  ,t)]Jr 
o<t} 


;=l,-,r 


2.4  Trajectory  Sensitivity  Function 

The  output  sensitivity  functions  defined  above  are  for  the 
systems  with  a  single  degree  of  freedom.  They  can  be  determined  by 
solving  the  high-order  ordinary  linear  differential  equation.  In 
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theory  the  same  approach  can  be  used  for  a  system  with  multiple 
degrees  of  freedom.  This  will  be,  however,  a  enormous  task  in 
practice.  Nevertheless,  the  governing  equation  can  be  conveniently 
written  in  terms  of  the  state  variables.  The  similar  approach  then 
can  be  adopted  in  the  determination  of  the  sensitivity  functions. 

Consider  a  system  govern  by  the  differential  equations 

x = F(x,u,t,a),  x(t0)=x$  (2-23) 


in  which  x  is  an  nXl  state  vector,  F  an  nXl  excitation  vector,  u  an 
mXl  input  vector,  and  a  an  rXl  parameter  vector.  Taking  the 
partial  derivatives  on  both  sides  of  equation  (2-23)  with  respect 
to  ctj  yields 


dx  _  dF  dx  |  dF  dx  , 

da  ■  dx  daj  da j  3ay.  ,=0 


7=1, -,r 

where  dF/dx  is  an  nXn  Jacobian  matrix, 


(2-24) 


fBFl  dF \y 

dx i  dxn 


dF 

dx 


(2-25) 


V-  ...  V. 

,  &1  &U) 

Again,  interchanging  the  derivatives  with  respect  to  the  time  t  and 
parameter  ot j ,  and  evaluating  the  above  equation  at  the  nominal  value 
a0,  we  arrive  at 
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(2-26) 


3a , 


7  =  1,  -  ,r 


where  £j  is  called  the  trajectory  sensitivity  function  with  respect 
to  the  j-th  component  of  the  parameter  vector  a,  and  is  defined  by 


(2-27) 

j 

Equation  (2-26)  generally  is  a  set  of  linear  time  variant 
differential  equations.  Since  coefficients  in  the  equation  contain 
nominal  solutions  x(at°,t),  the  numerical  computation  for  this 
problem  would  be  quite  expensive.  The  situation  will  be  very 
different  if  the  system  is  linear,  and  represented  by 

x=Ax  +  Bu;  x = xQ  (2-28) 


in  which  A  and  B  are  matrices  containing  a.  The  corresponding 
sensitivity  equations  appear  to  be  linear  time-invariant,  i.e., 


^=A^+^LoJC(M);  ^(0)  =  0 

7=1,  ~,r 


(2-29) 


2 . 5  Performance-Index  Sensitivity 

A  linear  quadratic  regulator  is  an  optimization  problem  of 
determining  the  optimal  control  force  u(t)  by  minimizing  the 
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performance  index 


J=G  +  f'jLdt 

where 


(2-30) 


G=x  T(t)Sx(tf)  (2-31) 

and 

m  =xT(t)Q(t)x(t)  +  u  T(t)R(t)  u(t)  (2-32) 

The  constraint  for  the  performance  index  is  given  in  equation  (2- 
28)  .  The  sensitivity  for  the  performance  index  is  a  vector  defined 
by 


T) 


—  I 

da 


(2-33) 


Taking  derivatives  of  equation  (2-30)  and  evaluating  it  at  the 
nominal  value  a0  yield 


(2-34) 


in  which  $  is  the  trajectory  sensitivity  vector  obtained  from 
equation  (2-26) . 

Methods  of  determine  other  kinds  of  sensitivity  functions, 
such  as  eigenvalue  and  eigenvector  sensitivity  functions,  can  be 
found  in  References  [3,4,5,15]. 
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3.  INTEGRATED  STRUCTURAL  CONTROL  AND  OPTIMIZATION 


Structural  design  is  of  an  exercise  to  aim  the  fulfillment  of 
safety  criteria  and  economical  expenses.  For  a  structural  system 
under  static  loadings,  it  is  a  common  practice  to  employ  the 
optimization  technique  to  minimize  certain  cost  function  of  the 
system,  such  as  the  strength  or  stiffness,  by  varying  system 
parameters  in  a  prescribed  range.  This  problem  can  be  stated  as  a 
linear  or  nonlinear  mathematical  program,  that  minimizes  the  cost 
function  J,(a) ,  where  a  is  the  vector  of  system  parameters,  with 
constraints,  ho  <  h(a)  <  ht.  Sensitivity  derivatives,  the 
important  ingredients  in  the  optimization,  can  be  easily  determined 
for  such  a  case.  When  a  structural  system  is  excited  by  dynamic 
loadings,  however,  the  fulfillment  of  the  safety  criteria  and 
economical  expenses  is  much  more  difficult  to  achieve  solely 
through  the  optimization,  since  either  the  cost  function  J,  or  the 
constraints  may  be  functions  of  time  t.  Attempt  has  been  made  to 
change  the  system  characteristics,  such  as  the  eigenvalues,  by 
varying  damping  ratio  and  stiffness  coefficients,  such  that  the 
dynamic  response  can  be  controlled  to  meet  some  specific 
requirements  [10,12].  This  effort  is  also  known  to  be  passive 
control,  a  limited  approach  especially  when  control  of  overall 
performance  of  the  structural  system  is  desired.  The  reason  for 
the  ineffectiveness  of  this  control  is  that  it  only  changes  the 
amount  of  energy  dissipation  or  the  resonance  frequencies,  but  by 
no  means  alternates  the  effect  of  the  input.  In  contrast,  the 
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active  control  provides  a  system  with  external  energy,  capable  of 
canceling  out  parts  of  the  energy  due  to  the  external  excitation, 
keeping  its  effect  on  the  dynamic  response  minimum.  This  method 
has  been  proved  to  be  very  efficient.  However,  experience  has 
shown  that  the  structural  control,  in  most  cases,  was  carried  out 
independently  from  the  structural  optimization.  Possible 
improvement  in  terms  of  combination  of  the  control  and  optimization 
had  long  been  overlooked,  until  recent  papers  [6,7]  that  drew  much 
attention  in  the  structural  engineer  community.  These  papers 
provide  a  clear  evidence  that  changes  of  system  parameters  can 
affect  control  performance  significantly.  They  strongly  suggested 
that  control  and  optimization  should  be  applied  to  a  structural 
system  simultaneously  in  order  to  achieve  the  best  structural 
design. 

It  is  of  importance  to  note  that  the  modern  control  theory  was 
established  on  the  basis  of  the  optimization  scheme.  As  an 
example,  the  linear  quadratic  optimal  control  is  to  minimize  the 
performance  index 


Jc(a,u)  =  j*f(x  tQx  +  uTRu)dt  (3- 1) 

with  constraints 

x  =  Ax+Bu;  x(0)  =x°  (3“2) 

in  which  x  and  u  are  vectors  of  state  variables  and  active  control 
forces,  respectively,  and  Q  and  R  are  appropriate  weighting 
matrices  determined  according  to  relative  emphasis.  Matrix  A 
contains  system  parameters  o,  which  in  general  for  a  structural 
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system  is 


A(a)  = 


0 

-M~\a)K(a) 


I 

-3 r\a)C(a) 


(3-3) 


where  M,  C  and  K  are  the  mass  matrix,  the  damping  matrix  and  the 
stiffness  matrix,  respectively.  The  dependence  of  a  in  the 
performance  index  is  due  to  the  constraint  which  in  fact  is  the 
governing  equation  for  the  structural  system.  For  optimal  control, 
optimization  is  performed  on  the  control  variables  u  with 
predetermined  system  parameters  a.  In  other  words,  if  the 
performance  index  Je  and  the  control  forces  u  are  treated  as  the 
system  response  and  the  system  parameters,  respectively,  the 
control  itself  is  a  optimization  process.  However,  the  control 
variables  are  usually  not  determined  by  the  conventional 
optimization  approach.  Instead,  they  can  be  obtained  indirectly  by 
solving  a  set  of  nonlinear  ordinary  differential  equations,  called 
the  Riccati  equation  [11], 


P+A  TP+PA-PBRlBTP+  <2  =  0  (3-4) 

in  which  P  is  a  matrix  attributed  to  the  optimal  feed  back  control 
forces 


u  =  -R~1BtPx  (3-5) 

Since  the  Riccati  equation  contains  the  coefficients  of  the 
governing  equation  explicitly,  the  original  optimization  is 
simplified  to  a  ordinary  differential  equation  problem. 

Many  researchers  suggested  that  the  integrated  structural 
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control  and  optimization  be  pursued  by  two  separate  paths 
[8,9,13,14]:  the  optimal  control  and  the  parameter  optimization. 
In  this  two-path  optimization  approach,  structural  design  variable 
a  and  the  control  variables  u  are  classified  as  two  sets  of 
parameters  to  be  optimized.  The  corresponding  objective  functions 
are  J,(a)  and  Jc(a,u) ,  respectively.  The  purpose  of  the  integrated 
structural  control  and  optimization  is  to  minimize  the  total 
objective  function,  J(a,u)  =  Js(a)  +  Jc(a,u)  ,  with  respect  to  a  and 
u  under  the  conditions  that 


h04h(a) 

x  =  Ax+Bu;  x(0)=xQ 


(3-6) 


Instead,  the  two-path  optimization  is  performed  on  each  sub¬ 
objective  function  separately.  The  minimization  of  J(a,u)  with 
respect  to  a  only  involves  J,(o) ,  but  not  J(a,u) .  The  algorithm  for 
this  path  of  optimization  is  exactly  the  same  for  the  passive 
control.  The  design  variables  a  are  fixed  while  the  J(a,u)  is 
minimized  with  respect  to  u.  In  this  manner,  the  total  objective 
function  is  minimized  to 


=  ^J/a)  +  n^nJc(a**,u) 


(3-7) 


The  second  part  of  the  minimization  usually  is  achieved  by  solving 
the  Riccati  equation,  which  is  the  main  advantage  of  this  approach 
as  far  as  the  numerical  computation  is  concerned.  Since  the 
coupling  of  the  two  sets  of  parameters  is  disregarded  in  each  path 
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of  optimization,  a”  and  u”  are  not  true  optimal  values,  and 
J(a**,u**)  is  not  the  true  minimum  value.  It  can  be  proved  that 

=/(«>•)  (38) 

Thus,  the  two-path  optimization  only  leads  to  sub-optimal 
solutions. 

Salama,  et  al.  [16]  proposed  an  alternative  in  which  the 
coupling  of  the  two  sets  of  parameters  are  taken  into  account.  The 
strategy  of  this  optimization  can  be  described  by  the  following 
expression 


Consider  the  special  case  where  the  control  is  activated  in  the 
period  from  zero  to  infinity.  x(t=°o)  =  o  then  is  realized. 

Recognizing  the  fact  that  [11] 

=  m^f~[xTQx  +  uTRu]dt  =  x0TP(  a)x0  (3-10) 

and 

u*  =  -RrlB(a)TP(a)x  (3-1 1) 

where  matrix  P(a)  is  the  solution  for  the  Riccati  equation  (3-4) . 
We  can  rewrite  equation  (3-9)  as 

J(oc’,u-).I^rl[J1(a)+i0TJ>(<t)xJ  (3-12) 

The  right  hand  of  the  above  equation  now  is  an  implicit  function  of 
a.  With  the  imposition  of  constraints  on  a,  the  optimal  solution 
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can  be  determined  by  using  sensitivity  analysis.  This  algorithm 
requires  the  derivatives  of  matrix  P(o),  which  must  be  obtained 
from  the  Riccati  equation.  Since  the  Riccati  equation  only  can  be 
solved  numerically,  calculating  of  the  sensitivity  derivatives  with 
respect  to  a  can  be  cumbersome.  The  detailed  discussion  of  this 
method  can  be  found  in  [16]. 

The  methods  stated  above  have  one  thing  in  common:  the  Riccati 
equation  being  used  in  the  optimal  control  part.  Although  solving 
the  Riccati  equation  requires  sophisticated  techniques,  this  effort 
would  be  still  much  less  than  that  using  the  conventional 
optimization  approach.  The  value  of  these  methods  largely  depends 
upon  the  legitimacy  of  the  Riccati  equation  for  specific  problems. 
In  what  follows,  we  will  carefully  reexamine  the  conditions  for 
which  the  Riccati  equation  is  valid. 

We  start  with  the  governing  equation  for  a  structural  system, 

x  =  Ax+Bu+Hf,  x(0)=xo  (3-13) 

where  f(t)  is  a  vector  of  external  forces.  The  objective  of  the 
active  control  design  is  to  minimize  the  performance  index  in 
equation  (3-1) .  As  stated  early,  equation  (3-13)  is  of  the 
constraints  for  the  optimal  problem.  Now  we  combine  equations  (2- 
1)  and  (2-13)  and  form  a  Lagrange 

c£  =  tQx  +  utRu-  X  t(x-Ax -Bu-  Hf)]dt  (3’14) 

where  X  is  a  vector  of  Lagrangian  multiplier,  which  in  general  is 
time  dependent. 
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Taking  variation  of  equation  (3-14)  results  in 


6S£=  ft'{(XT+—) 5x+  —}dt+XT(0)6x(0)  -  (3-15) 

Jo  dx  du  11 

where  &  is  the  integrand  in  equation  (3-14)  .  The  necessary 
condition  for  ££  to  be  minimum  is  SX  =  0,  which  corresponds  to 

—  =0  (3-16) 

du 

X  +  —=0  (3-17) 

dx 

and 

1(9  =  0  (3-18) 

The  second  term  on  the  right  hand  side  of  equation  (3-15)  vanishes 
due  to  the  initial  condition.  Equations  (3-16)  and  (3-17) 
constitute  a  new  set  differential  equations  of  X, 

X  =  -AX-2Qx\  1(9  =  0  (3-19) 

For  a  closed-loop  control,  X(t)  is  taken  to  be 

X(t)  =  P(t)x(t)  (3-20) 

The  corresponding  optimal  control  force  can  be  determined  according 
to  equation  (3-16) , 

u(t )  =  -R-WXCt)  (3-21) 

Substituting  it  into  equation  (3-17) ,  we  obtain 
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(P+ATP+PA- PBR  ~1BtP+Q)x+ PHf  =  0 

(3-22) 

O 

II 

(3-23) 

For  the  special  case  where  f  =  0,  which  usually  happens  for 
parametric  excitations  or  structure-induced  excitations,  equation 
(3-22)  reduces  to  the  Riccati  equation  (3-4)  .  It  also  can  be 
proved  [11]  that  the  optimal  solution  for  the  case  in  which  f  is  a 
white  noise  excitation  has  the  same  form  of  equation  (3-21) ,  and 
corresponding  matrix  P  again  can  be  solved  from  the  Riccati 
equation  (3-4) .  Therefor,  the  necessary  condition  for  the  Riccati 
equation  to  be  valid  for  the  optimal  control  is  that  the  additive 
external  excitation  f  is  either  white  noise  or  zero.  Since  most 
excitations  applied  to  the  structural  systems  are  not  white-noise, 
the  optimal  control  cannot  be  determined  from  the  Riccati  equation. 
As  the  consequence,  theories  presented  early  for  the  integrated 
control  and  optimization  have  limited  applications  for  structural 
systems . 

Dealing  with  a  structural  system  subjected  to  a  general 
loading,  Hale  and  his  colleagues  [7]  derived  the  necessary 
conditions  for  optimal  control  and  optimal  system  parameters  by 
using  variation  principle.  The  objective  function  they  used  is 

J  =  f*'  [xtQx  +  u  tRu  +Js(a)-X  T(t)[x  -  A(a)x  -  B(a)u  -  H(a]f]]dt  (3-24) 
in  which  Jc(a)  is  the  cost  function  for  design  variable  a.  The 
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condition  for  the  minimum  J  can  be  obtained  by  taking  variation  on 
equation  (3-24) .  These  conditions  are 


x=Ax+Bu+Hf;  x(0)  =  0 

(3-25) 

2Qx+ATX  +  X=0 ;  1(^  =  0 

(3-26) 

2Ru+BTX=0 

(3-27) 

X  rC Vjix  +  VJBu  +  VaHf)  +  VB  W  =  0 

(3-28) 

It  is  of  importance  to  note  that  the  advantage  of  this  approach  is 
that  it  directs  the  solution  toward  to  a  single  optimal 
configuration  so  that  the  true  optimal  design  of  control  forces  and 
system  parameters  is  achieved.  Bearing  it  in  mind  that  matrices  A, 
B  and  H  are  functions  of  parameter  a,  the  above  equations  are 
coupled  in  a  and  u.  It  is  almost  impossible  to  solve  them 
simultaneously  for  analytical  solutions.  Therefore,  numerical 
computations  are  necessary  to  search  the  optimal  solutions. 
Several  computational  technique  are  available  to  solve  equations 
(3-25)  through  (3-28)  [7].  They  appeared  to  be  very  time 
consuming.  In  addition,  they  have  not  been  fully  evaluated  in 
terms  of  the  rate  of  convergence. 
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4.  OPTIMIZATION  BASED  INSTANT  OPTIMAL  CONTROL 


As  mentioned  before,  the  theory  of  the  optimization  is  well 
developed  along  with  some  standard  optimization  algorithms. 
However,  increasing  applications  of  large  scale  structures  in  the 
modern  time  demand  further  developments  toward  to  the  efficiency  of 
optimization  techniques.  Recent  advances  shows  that  there  have 
been  great  achievements  in  such  a  mission.  One  of  important 
developments  in  optimization  algorithm  attributes  to  a  recent  work 
by  Venkayya  [20],  who  developed  a  so-called  compound  scaling 
algorithm,  which  has  been  proved  to  be  very  efficient  for  large 
scale  structures.  The  objective  of  current  research  is  to  extend 
this  new  method  to  optimal  control  problem,  and  possibly  to  the 
integrated  control  and  optimization  for  active  structures.  In 
what  follows,  we  will  propose  a  new  control  algorithm  in  which  the 
optimal  control  can  be  completely  treated  as  static  parameter 
optimization. 

The  conventional  optimal  control  relying  on  the  Riccati 
equation  requires  that  the  external  excitations  either  be  zero  or 
white  noise.  The  nature  of  the  excitation,  in  reality,  however,  is 
random  and  usually  cannot  be  predicted  exactly.  Hence,  the  control 
design  must  depend  upon  the  instant  information  of  the  excitation, 
and  should  be  implemented  instantly.  To  this  end,  we  formulate  an 
instant  performance  index 

J(t)  =  *  T(t)Q(t)x(t)  +  u  T(t)Bx(t)  (4'1) 
Our  purpose  is  to  minimize  J(t)  at  every  time  instant  with 
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constraints 


x(t)=Ax(t)+Bu(t)+Hf(t);  x(0)  =x0  (4-2) 

For  a  linear  system,  it  can  be  proved  that  over  a  small  time 
interval  At,  the  state  vector  can  be  written  as 

x(f)  =  Tdit-M)  +  -^[Bu(t)+Hf(t)]  (4-3) 

mm 

where  T  is  a  transformation  matrix  independent  of  time  t,  and  d(t- 
At)  is  a  vector  containing  previous  information  of  x(t-At)  ,  u(t-At) 
and  f(t-At).  Equation  (4-3)  is  equivalent  to  equation  (4-2),  and 
again  is  a  constraint  for  the  instant  performance  index.  For  the 
closed-loop  control,  the  control  force  is  directly  related  to  the 
state  vector  x(t) ,  namely; 

u(t)=P(t)x(t)  (4-4) 

in  which  P(t)  is  a  control  gain  matrix.  Specific  requirement  can 
be  imposed  to  this  control  gain  by 

W(f)]-9t0(')  =  0  (4-5) 

is  a  operator  which  may  represent  eigenvalues,  eigenvectors, 
determinant  or  other  properties  of  matrix  P(t).  Equation  (4-5) 
along  with  equation  (4-4)  are  additional  constraints  for  the 
instant  performance  index. 

Equations  (4-1)  through  (4-5)  constitute  a  regular  parameter 
optimization  problem.  They  can  be  cast  into  a  Lagrange 
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(4-6) 


a =7(o -21  (fix/®  -  s  tfimt) -  z  |i  kd)zt(t) 

1  10 

where  Xj(t)  ,  7j(t)  and  nk( t)  are  Lagrangian  multipliers  corresponding 
to  the  constraints 

X(t)  =  *(f)- (4-7) 


Y(t)  =  u(t)-mx(t)z  o 


(4-8) 


and 


Z(f)  =  M-3t0(^0 


(4-9) 


In  equation  (4-6) ,  active  variables  participating  the  optimization 
are  x(t) ,  u(t)  and  P(t),  which  can  be  cast  into  a  qXl  vector  a. 
The  stationary  conditions  for  the  Lagrange  are 


dg(t)  _  dJ(t) 
3a p  dap 


-SX,(0 

1 


-sYj(0 

i 


dY/t) 

da 

p 


n 


-S  »k(t) 


dZfi) 

dap 


=  0; 


(4-10) 


p  =  \,2,-,q 

Equation  (43)  then  can  be  written  as  the  Kuhn-Tucker  conditions 


2Je..A  .+Hh..\ .+ Ejp.li.  =  1; 

i y  j  1 1  j  6ij^j  5 

1  1  X  J 


i=l,2,-,q 


(4-11) 


where 
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(4-12) 


dXj 

da . 

e.  =  — - 

*  a/ 


aa. 


ii. 

a/ 

3ay. 


(4-13) 


dz. 

da. 

f?.- - J- 

JU_ 

da . 
j 


(4-14) 


are  ratios  of  sensitivity  derivatives,  which  can  be  determined 
using  the  methods  in  Chapter  2. 

It  is  worthy  to  note  that,  at  any  time  instant,  equation  (4-7) 
can  be  viewed  as  the  conditions  for  static  equilibrium.  The 
formulations  presented  in  equations  (4-6)  through  (4-9)  constitute 
a  standard  parameter  optimization  problem  at  any  time  instant  t, 
provided  that  the  optimal  solutions  x(t-At)  and  u(t-At)  are 
available.  This  requires  that  the  optimal  solution  x(t)  and  u(t) 
be  obtained  within  a  time  period  At  in  order  to  progressively  solve 
the  problem  in  a  full  time  scale.  The  on-line  computation  must 
engage  if  better  accuracy  is  required.  Therefore,  the 
computational  efficiency  is  essential  to  the  proposed  scheme. 
Among  the  existing  numerical  techniques,  the  compound  scaling 
algorithm  [20]  is  the  most  feasible  one  for  the  current 
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optimization  based  instant  control.  This  method  was  developed  by 
successfully  applying  the  scaling  and  resizing  techniques  in  the 
optimization  procedure.  The  remarkable  feature  of  this  approach  is 
that,  instead  of  searching  the  optimum  point  by  point,  like  the 
conventional  ones,  in  the  design  space,  the  search  sweeps  the 
design  space  and  approaches  the  optimum  in  a  fast  rate  [20]. 

The  optimization  based  instant  optimal  control  has  explicit 
superiorities  over  the  conventional  optimal  control.  First,  it  can 
be  applied  to  structural  systems  under  general  dynamic  loadings. 
It  is  especially  useful  when  the  loading  is  random.  Secondarily, 
The  performance  index  does  not  contain  integrations  or  derivatives, 
and  optimal  control  is  designed  for  time  instant  t.  Therefore,  the 
problem  is  suitable  to  be  treated  as  regular  static  structural 
optimization.  Finally,  the  optimal  control  solutions  satisfy  the 
Kuhn-Tucker  conditions.  The  compound  scaling  algorithm  can  be 
used;  so  that  numerical  computation  is  expected  to  be  very 
efficient. 


30 


5.  CONCLUDING  REMARKS 


In  this  report,  we  have  presented  several  methods  to 
indirectly  determine  sensitivity  measures  related  to  structural 
controls  in  the  time  domain.  Upon  reviewing  the  existing 
techniques  in  the  area  of  integrated  control  and  optimization,  we 
conclude  that  approaches  that  relies  on  the  Riccati  equation  are 
limited  to  the  structures  subjected  to  zero  additive  excitation  or 
white  noise  random  excitations.  True  optimal  solutions  should  be 
obtained  by  treating  design  variables  and  control  variables  equally 
as  the  system  parameters  to  be  optimized.  The  sensitivity  analysis 
obviously  should  play  the  key  role  in  this  approach.  Future 
research  should  concentrate  on  the  development  of  efficient 
computational  schemes  that  can  implement  the  sensitivity  analysis 
in  the  integrated  control  and  optimization  systematically.  Also, 
we  proposed  a  new  optimal  control  algorithm,  capable  of  converting 
a  optimal  control  problem  to  a  static  parameter  optimization 
problem.  This  new  algorithm  provides  us  an  opportunity  to  take  the 
advantage  of  recent  advances  in  the  optimization  area.  This  method 
has  potential  to  be  extended  to  active  structures  by  including 
system  parameters  in  the  active  parameters  in  optimization.  Better 
performance  of  the  active  structures  can  be  achieved  due  to  the 
fact  that  the  structural  properties  and  control  variables  are 
optimized  at  every  time  instant.  Since  this  report  is  prepared  for 
our  preliminary  research,  many  details  have  not  been  touched 
thoroughly.  Prospective  results  will  be  developed  in  the 
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